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Abstract. Being chosen as a differential operator of a special form, metric 77 
operator becomes unitary equivalent to a one-dimensional Hermitian Hamiltonian 
with a natural supersymmetric structure. We show that fixing the superpartner of 
this Hamiltonian permits to determine both the metric operator and corresponding 
non-Hermitian Hamiltonian. Moreover, under an additional restriction on the non- 
Hcrmitian Hamiltonian, it becomes a superpartner of another Hermitian Hamiltonian. 
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1. Introduction 



■ After a seminal paper by Scholtz Geyer and Hahne [1] , metric operator 77 plays a crucial 

role in the pseudo-Hermitian (quasi-Hermitian) quantum mechanics. Just this operator 
{Nj ■ is used while one redefines the inner product to bring initially non-Hermitian operator 

to its Hermitian form (for details see review papers [2], [3]). Usually (see e.g. [3]) this 
is an invertible, Hermitian, positive definite and bounded operator such that 

7]H = H^T]. (1) 

Here H is a given non-Hermitian operator (Hamiltonian) and W is its Hermitian adjoint. 
If H has a discrete spectrum and admits a complete set of biorthonormal eigenvectors 
then the metric operator may be presented as [4] 

7] = 00^ . 

Note that this form is typical (see e.g. [5, 6]) when a Hermitian Hamiltonian is 
presented in a factorized form in supersymmetric quantum mechanics (SUSY QM). 
Its superpartner r] is obtained by interchanging O and 

Vo = 0^0. 

The only difference between the SUSY intertwiner, that we will denote L, and the 
operator O is that L is a differential operator so that the metric LU becomes a 
differential operator. It is clear that LL^ is Hermitian. It is both positive definite 
and invertible if it has an empty kernel in the corresponding Hilbert space. This is 
possible if both L and have empty kernels. In SUSY QM this corresponds to a 
broken supersymmetry [6]. 
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Thus the only obstacle to use the technique of SUSY QM for studying properties 
of differential metric operators is that these operators, being differential operators, 
are usually unbounded. Nevertheless, for instance Fityo [7] used a first order 
differential operators O for constructing a new class of non-Hermitian Hamiltonians 
with real spectra without any discussion that the corresponding metric operator becomes 
unbounded. He also mentioned that H and W are superpartners of the second order 
supersymmetry. Nevertheless, in contrast to the usual second order supersymmetry 
[5], here an additional restriction is imposed on the intertwining operator, it should 
be Hermitian. Consequences of this restriction are not analyzed in [7]. Our analysis 
shows that corresponding first order operator L (actually its complex conjugate form 
L*) intertwines H with a Hermitian operator that we will denote h . 

Any unbounded operator has a domain of definition which is a subset of the 
corresponding Hilbert space. Therefore while redefining the inner product in the spirit 
of paper [8] 

«^i|^2», = <^ihl^2> , ^1,^2 en (2) 
where H is the initial Hilbert space, one has to replace the condition V^i,^ G % by 
another one V0i, ^2 G D v C H where D v is the domain of definition of i] operator. From 
here it follows that there exist in % vectors tp ^ D v which cannot be mapped into the 
new Hilbert space H v defined with the help of the inner product ((V'iIV^))??- This is the 
price for using unbounded rj operators. From the physical viewpoint, it remains to hope 
that maybe states described by such vectors cannot be realized in practice. Moreover, 
this is not an obstacle for finding Hermitian Hamiltonian h related to the given non- 
Hermitian if by a similarity transformation [9], [10]. In particular, in [10] scattering 
matrix and cross section for h are calculated and their unusual properties are discussed. 

2. Second order differential rj operator 

Let we are given a non-Hermitian Hamiltonian 

d 2 

H = -— + V(x), x>0 (3) 
dx z 

with the domain 

D H = {<^GL 2 (0,oo) : (p"(x) — V(x)tp(x) G L 2 (0, oo) , [ip\x) + w(x)<p(x)] x=0 = 0} (4) 

where V(x) and w(x) are complex valued functions. We find adjoint operator W as 
usual using the inner product un the space L 2 (0, oo), (Wipi\(p 2 ) = (<pi\H(p2), y?2 £ Dh, 
<fi G D H i thus obtaining 

D Ht = {y? G L 2 (0, oo) : <p"(x) - V*(x)<p(x) G L 2 (0, oo) , [<p'(x) + w*(x)ip(x)] x=0 = 0} 

& = -^ + V*(x) = H*. (5) 
dx z 

For simplicity we will assume that V(x) is a scattering potential satisfying the condition 

(l + a;)|V(a;)|cZa;<oo. (6) 
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In particular, we will assume that the function \ V(x) | is bounded below and tends to zero 
faster than any finite power of 1/x when x — > oo so that the operator of multiplication 
by the function V(x), as an operator acting in L 2 (0,oo), becomes bounded. For this 
reason, we can simplify the domains Dh and D H \ 

D H = { V EL 2 (0,oo): <p"(x) GL 2 (0,oo), [<p'(x) + w(x)ip(x)] x=0 = 0} , (7) 
D H t = {y? G L 2 (0, oo) : <p"(x) G L 2 (0, oo) , [<p'(x) + w* (x)<p(x)] x=0 = 0} 

As was discussed in Introduction, to find the metric operator for the given non- 
Hermitian Hamiltonian H, one has to find a Hermitian positive definite and invertible 
operator 77 satisfying equation (1). If 77 is bounded its domain is the whole Hilbert space 
and no problems occur to act both by the left and by the right hand sides of (1) on 
functions belonging to Dh- Unfortunately, this is not our case since we want to consider 
unbounded 77 having its own domain in L 2 (0, 00). We find resonable to assume that the 
domain of 77 coincides with that of H, D v = Dh (4) (or (7)). As we show below, if 
operator 77 is chosen to be a second order differential operator of a special form defined 
on this domain, this assumption is justified since such 77 is selfadjoint. Nevertheless, 
even in such a case one cannot apply (1) to any ip G D H . Indeed if ip G D v — D H then 
from the left hand side of (1) it follows that Hip G D H . Thus (1) has a sense on a subset 
of Dh such that the range of H is contained in the domain of H . It is possible to show 
that this subset is dense in L 2 (0, 00) but we do not dwell on its proof. 

Assume that n is defined with the help of a differential expression. Then from the 
condition D v = D H it follows that n may be a second order differential operator for 
which we assume the form 

77 = LL t (8) 

where 

L = —d/dx + w*(x) , L j = d/dx + w(x) . (9) 

so that 

nip(x) = -ip"(x) + (w* - w)ip'(x) + [ww* - w']ip(x) , ip(x) G D rj = D H ■ (10) 

We note that operators L and are mutually adjoint with respect to the inner product 
in L 2 (0, 00). To see that consider for instance 

rco poo 

(^\l^ 2 )= rM + ™4>2)dx = [rM%>+ (-4> > 1 + w*)*i> 2 dx = (Li> 1 \4> 2 ). 

Jo Jo 
The integrated term here vanishes at infinity since the functions ip\{x) and ip 2 {x) are 
smooth enough and belong to L 2 (0,oo) and therefore they tend to as x — > 00. At 
x — it vanishes since ipi G Dl and the domain of L coincides with the range of 
which according to (7) reads 

D L = {^p G L 2 (0, 00) : ip'(x) G L 2 (0, 00), ^(0) = 0} . (11) 

It is easy to check that the operator 77 (10) is selfadjoint with respect to the inner 
product in L 2 (0, 00). Indeed, as usual assuming that ip\ G D v and integrating by parts 
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twice the term with the second derivative and once the term with the first derivative 
yields 

POD 

Jo 

= [ ^2^1 - $2^1 + ^2^1 ( W * ~ W ) }7=0 

POD 

+ / [ -^2*" + (w- w*)i)* 2 ' - (w*' - w')^l + (ww* - w')ip* 2 ] i) X dx 
Jo 

POD 

= I [—lp2 + (w*—w)lp2 + ( ww *~ w ')*l J 2}**Pldx 

Jo 

= (v^tpi) ■ 

To justify the last equality we consider the integrated term at x = 

[^Vl - + V>2*V>1 {W* - W) } x=0 = [^2*Vl - V>2*(-W>l) + {W* - W) } X=Q 

= iMo)[^' + ^xU> 

= 0. 

The first line here follows from the property that ip\ G D v = Dh given in (7) 
and in the last line we used ip2 G D q = Dh- In Section 7 we consider examples 
illustrating selfadjointness of i] when it is a second order differential operator with 
constant coefficients (w'(x) = 0) and with variable coefficients (w'(x) ^ 0). 

Furthermore, if we impose the condition that the operator has the empty kernel 
in D v then i] is positive definite. Imposing additionally that the operator L has the 
empty kernel in Dl (11), we get an invertible operator rj. With these assumptions, we 
have operator r) suitable for constructing operator p = r^l 2 (see below). 

Another property that follows from the condition D v = Dh (4) (or (7)) is that the 
differential equation 

r]^ k (x) = \(k)* k (x) , k > (12) 

should have only one singular point which is x = 00. This means that the function w(x) 
should be regular for all x G (0, 00) and therefore we can put 

w (x) = dlo ^(x)} ^ u ( x ^ VxG(0,oo). (13) 
dx 

Since D v = Dh (7) we have to supply equation (12) with the following boundary 
condition at x = 

[%(x)+w(x)V k (x)} x=0 = 0. (14) 

For simplicity we will consider the case when operator rj has no bound states. Therefore 
we will impose on the functions ^(x) an asymptotic condition at x — > 00 so that they 
describe scattering states of the operator 77. 



3. Unitary equivalence between rj operator and a Hamiltonian 

Let us put 

u{x) = p{x)e iw{x) (15) 

where p(x) and ou(x) are real valued functions. Then after a unitary transformation 

^ k = U^ k , n = UWA ] (16) 

with the operator U being a multiplication operator on function (actually it is simply a 
phase factor) 

«=(ff--"' (if = .* (17) 

equation (12) reduces to a Schrodinger-like equation 

H* k (x) = \(k)* k (x), U = -^+V(x) (18) 

where the potential V(x) is defined in terms of a real valued function (superpotential) 
W(x) as follows 

V(x) = W 2 (x) - W'(x) (19) 
and the function W(x) is expressed in terms of the modulus of the function (15) 

W(x) = \ [log(««*) ]' = [logp]' = ^ = W(s) . (20) 

Boundary condition for equation (18) follows from that for equation (12) 

[*' k (x)+W(x)* k (x)] x=o = (21) 
and from (4) we get 

D H = {* G L 2 (0, oo) : V"(x)-V(x)*(x) G L 2 (0, oo) , [&(x)+W(x)*(x)] x=0 = 0} .(22) 

Potential (19) has the structure typical for the supersymmetric quantum mechanics [5]. 
Its SUSY partner Vo has the form 

V (x) =W 2 (x) + W(x) = ^ . 

p{x) 

and for the potential V(x) we obtain a formula typical for the Darboux transformed 
potential (see e.g. [5]) 

V(a;) = V (x) - 2W'(x) = V (x) - 2 [logp(ar) ]" . (23) 

Moreover, the function p(x) is a solution to the Schrodinger equation with the 
Hamiltonian 

Ho-^+Vo(.) (24) 
corresponding to zero eigenvalue 

H p(x) = . (25) 
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Function p(x) is known as a transformation function for the Hamiltonian H (see e.g. 
[5]). If solutions to the Schrodinger differential equation with the potential ~Vo(x) are 
known 

H ^ ) (x) = E^\x) (26) 

then solutions ^e(x) to the same equation with the potential V(x) may be obtained 
with the help of the transformation operator 

L p = -d/dx + W(x) = -d/dx + p'{x)/p{x) (27) 

as 

* E (x) = L p *%\x) = -[*%\x)]' + W(x)* ( §\x). 
Hence, Hamiltonians H and H are intertwined by the operators L p and L p 

L p U = iiL p , LjH = H Lj, (28) 

where 

L\ = d/dx + W(ar) . (29) 

From the first relation (28) follows that L p transforms eigenfunctions &^\x) of the 
Hamiltonian Ho (18) to eigenfunctions *&k(x) of the Hamiltonian H 

* k (x) = \^(k)L p ^\x). 

From the second relation (28) follows that L p transforms eigenfunctions *&k( x ) °f the 
Hamiltonian H (18) to eigenfunctions &f\x) of the Hamiltonian H 

*£\x) = X-^L^ix) = X-^WWix) + W(x)9 k (x)] , (30) 

H *f\x) = \(k)^\x). (31) 

Comparing equation (30) with (21), we find boundary conditions for the 
Schrodinger equation with the Hamiltonian H , 

*1 0) (0)=0. (32) 

It is easy to check that function 

u(x) = ix) (33) 

satisfies both the equation 

Hu(x) = (34) 

and the boundary condition (21). 

Operator H is unitary equivalent to 77. Therefore to fulfill the assumption that rj 
has a continuous spectrum without bound states, we have to assume that the operator 
H also has purely continuous spectrum and the Darboux transformation with the 
transformation function u(x) (33) over the Hamiltonian H does not create bound 
state. Furthermore, we have to assume also that the Darboux transformation from the 
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Hamiltonian H to the Hamiltonian H realized with the help of the transformation 
function p(x) does not create bound state either. Assuming additionally that the 
potential Vq(x) is scattering, we have to choose solution to equation (34) such that 

p^e dx , x^oo, d<0. (35) 

The choice d < 0, which we accept for what follows, guaranties an increasing asymptotic 
behavior of the function u(x) (33). Therefore this function is not square integrable and 
X(k) =0 does not belong to the spectrum of H. In this case operators rj, H and H are 
isospectral and corresponding supersymmetry is broken. We would like to emphasize 
that no any restriction on the phase ou(x) of the function u(x) (15) is imposed. 



4. SUSY partner for r] operator 

Operators L and D (9) play a crucial role in our approach since they define both 
rj operator and, as we show below, the Hamiltonians H and H^. Operators (9) are 
uniquely determined when the function u(x) is fixed (see equations (9) and (13)). 
Therefore properties of 77 operator and Hamiltonians H and depend on properties of 
the function u(x). Note that the modulus of u(x) is defined by the Hamiltonian H as 
a solution to equation (25). Therefore both the form of the operator 77 and properties 
of the Hamiltonian H depend on properties of the Hamiltonian Ho- 

Equation (12) may be considered as equation (18) unitary transformed with the help 
of operator U^ 1 (17) and therefore, as we show below, it inherits all supersymmetric 
properties of equation (18). 

Applying operator V to both sides of equation (12) yields, 

(L^L)L^ k (x) = X(k)L^ k (x) . (36) 
Above discussed condition X(k) 7^ implies 

*<®(x) = lS* k (x)?0. (37) 
From here one deduces that 

Vo ¥ k °\x) = X(k)¥ k °\x) , Vo = L^L. (38) 
Operators 77 (8) and t/q (38) are intertwined by V< and L, 

Vo tf = iSr) , L Vo = V L . (39) 

Thus operator L' maps solutions of Eq. (12) (eigenfunctions of rj) to solutions of Eq. 
(38) (eigenfunctions of 770 which is a SUSY partner of 77). Quite similarly, operator L 
realizes an inverse mapping. The transformation function for each mapping is either 
that which is annihilated by the operator L or that which is annihilated by its adjoint 
Z/t. The mapping preserving the ^-function normalization of the functions and 
is given by 

tf fc (a;) = X- l ' 2 (k) mf{x) , ^f\x) = X-^ik) lS* k (x) . (40) 
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Note that intertwining relations (39) are nothing but the associativity of operator 
multiplication 

(l)L)l) = l)(Ll)) L{pL) = (Ll))L . (41) 

We can resolve intertwining relations (1) in a similar way. For that we need not only 
operators L and V* but also their complex conjugate form 

L* = —d/dx + w(x) , (L*y = d/dx + w*(x) . (42) 

We assume here that the operation of the complex conjugation commutes with the 
operation of the Hermitian conjugation. 
Let us put 

H = L*tf + a, H^ = L(L*y + a* (43) 

where a is a complex constant. Then under an additional assumption 

L^L* + a = (L*) f L + a* (44) 

we reduce equation (1) to the identity. We thus expressed Hamiltonians H and W 
in terms of the function w(x) and, taking into account formula (13), in terms of the 
function u(x). We would like to emphasize that if the Hamiltonian H is fixed then the 
absolute value p{x) of the function u(x) is determined from equation (25) but its phase 
u(x) still remains arbitrary. Below we show that this arbitrariness may be fixed with 
the help of condition (44). 



5. Fixing phase oo(x) with the help of a Hermitian Hamiltonian ho 

Note that the right hand side of equation (44) is Hermitian conjugate with respect to 
its left hand side. This means that this equation becomes identity if operator 

h = L^L* + a = — ^- + w 2 (x) + w\x) + a = + + a ( 45 ) 

dx z dx l u[x) 

is Hermitian, h§ = h\. The necessary condition for that is the reality of the function 



v (x) = + a = ^j4- [uj\x)f + P + t L"{x) + 2^- u'(x) + 7 ) (46) 

u(x) p[x) V PV X ) 

where we put 



a = f3 + ry (47) 
and used equation (15). From here we find the equation for the function ou(x) 

W "(i) + 2^^)+ 7 = (48) 
and the potential 

v (x) = ^-lu'(x)f + (3 = v*(x) (49) 
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defining Hermitian Hamiltonian 

d 2 

ho = + v (x) = hl = h* . (50) 

From formula (45) we extract two important consequences. First, comparing it with 
(43), we conclude that the operators h and H are intertwined by operator L* (42), 

L*h = HL* . (51) 

Second, function u(x) is an eigenfunction of ho, 

hou(x) = au{x) . (52) 

Now we can formulate two approaches for finding a pair of operators H and r\, 
satisfying equation (1). Both approaches are based on the existence of an exactly 
solvable Hermitian Hamiltonian. 

In the first approach it is assumed that we know solutions to equation (52) with the 
given potential vq(x) both as a differential equation and as a spectral problem on the 
space of smooth enough functions from L 2 (0, oo). To be consistent with the previous 
assumptions imposed on the Hamiltonians H and H (see, e.g., equation (6)), here we 
will assume that v (x) is a scattering potential and ho has a purely continuous spectrum 

h ipk{x) = k 2 ip k {x) , k>0. 

Then taking a nodeless complex valued solution to equation (52) (parameter a may 
be both real and complex), we construct, with the help of equations (42) and (13), 
transformation operator L*. Moreover, from equation (51) we conclude that the 
Hamiltonian H (3) is a SUSY partner of ho (45) and therefore 

V(x) = v (x) -2[logu(x)\" . (53) 

Operating with the operator L* (42) on the eigenfunctions of the Hamiltonian ho, 
we obtain eigenfunctions of the Hamiltonian H . Operator 77 is found from equation 
(8). To find Hermitian operator h equivalent to non-Hermitian H, we have to solve 
eigenvalue equation (12) for the operator 77. For this purpose eigenvalue equations for 
the Hamiltonians H (18) or H (24) may be useful. 

Using intertwining relation (51) and properties ho — h\ — h* Q and W = H* (see 
equations (49) and (5)) yields 

h (L*y = (L*yH. 

This relation means that operators ho and H are SUSY partners and the operator (L*y 
(42) transforms eigenfunctions <Pk(x) of the Hamiltonian H to eigenfunctions ipk( x ) of 
the Hamiltonian h [10] 

Mx) = (k 2 - a) 1/2 [if' k (x) + w(x)ip k (x)\ . 

Comparing this equation with the boundary condition (7) for the eigenfunctions <fk{x) 
of the Hamiltonian H , we find the boundary condition for the eigenfunctions ipk{x) of 
the Hamiltonian h 
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As was mentioned in Introduction, we are interested to have a broken supersymmetry. 
Therefore we will consider the case when H and ho are isospectral. Hamiltonian H 
is assumed to have only continuous spectrum, hence, ho also should have no bound 
states and the transformation function should coincide with the Jost solution for the 
Schrodinger equation with the Hamiltonian ho- Then taking into account equation (52) 
and condition (35) yields the asymptotic behavior of the transformation function u(x), 

u{x) ->• e (d+ib)x , d<0, a = -{d + ib) 2 , (54) 

where b is an arbitrary constant. From here and equation (47) one finds the relation 
between constants d, b and j3, 7 (see (46)) 

(3 = b 2 - d 2 , 7 = -Idb . (55) 

This approach is realized in [10]. 

In the second approach the starting point is the spectral problem (31), (32) with 
the given potential ~Vq(x). Absolute value of the transformation function is found from 
equation (25) and its phase is fixed by equation (48). Potential V(x) follows from 
equation (23) and r\ operator is fixed by equations (16) and (17). Once both the modulus 
and phase of the transformation function (15) are fixed, we reconstruct Hamiltonian H 
(43) and if necessary potential v (49) and Hamiltonian h (50). We illustrate this 
approach in Section 7 by a simple example. 



6. Equivalent Hermitian Hamiltonian h 

Once eigenf unctions ^ k (x) and eigenvalues X(k) of the operator rj are found, one can 
write down its spectral representation 



V = / dk\(k)\* k )(9 k \. (56) 
Jo 

It has a unique Hermitian, positive definite and invertible square root 

p = V /2 =P f , P>0. (57) 

POO 

p = / dkX 1/2 (k)\ty k )(fy k \ . (58) 
Jo 

Hence, from equation (1) one finds the Hermitian operator h equivalent to H, 

h = pHp- 1 = p~ l H^p = (59) 

which in our case also has a purely continuous real and non-negative spectrum. Its 
eigenfunctions 

h$ k = k 2 $ k , k>0, 
are obtained by applying the operator p to the eigenfunctions of H [10] 

$ fc = (k 2 - a)-^ 2 PVk = (k 2 - a)"V^>fc • (60) 
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The factor (k 2 — a) x l 2 is introduced to guarantee both the normalization of these 
functions, 

= 5(k-k'), (61) 

and their completeness 

dk\$ k )($ k \ = 1. (62) 

From Eq. (60) it follows that the eigenfunctions (p k of h and $ fe of h are related by an 
isometric operator U [10] 

= 1% (63) 

where 

U = pL*(/i - a)' 1 = L(L t L)- 1 / 2 = L % 1/2 = (U^ 1 . (64) 

Using the spectral representation of the operator h, one can express h in terms of 
p, L*, (L*Y and the resolvent of h [10] 

h = [a(h - a)' 1 - a*(h - a*)' 1 ] (L*) f p . (65) 

Note that for any a = d + ib with d < the point a = —a 2 does not belong to the 
spectrum of ho and the operator (65) is well defined. When a is real a = a r = a* this 
expression becomes undetermined. This indeterminacy may be resolved using the usual 
l'Hospitale rule, which yields 

h = pL*h {h -a r )- 2 (L*y P . (66) 

From here we extract an important consequence. The real character of a r implies that 
d = and, hence, the point E = k 2 = a r = b 2 belongs to the continuous spectrum of H 
and corresponds to the spectral singularity [10]. This point belongs to the spectrum of 
ho also. For this reason the resolvent of h diverges at k 2 = a r = b 2 and operator 
(66) becomes undefined in the Hilbert space. This means that the non-Hermitian 
Hamiltonian H with a spectral singularity does not have Hermitian counterpart. 

Another possibility would be b = 0, d ^ 0. In this case a = —d 2 , operator H 
becomes Hermitian and h becomes unitary equivalent to H. 

Operator U (64) mapping eigenfunctions of ho to that of h (63) may be written in 
terms of the eigenfunctions \l/^ ' 1 of equation (38) 



poo 

U = L* dk\- 1/2 (k) |*i 0) )(* 
Jo 



(0), 
k I 



(67) 



Although equation (60) (or equations (63), (67)) formally solves the problem of finding 
the eigenfunctions of h, it contains the non-local operator p (or (L^L) -1 / 2 ) and, therefore, 
in general, no explicit expression for exists. 
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7. Examples 

As discussed in Section 5, there exist two approaches for finding a pair of operators 
H and r\ satisfying equation (1) and below we give two corresponding examples. First 
we exemplify the second approach since it leads to the simplest form of r] being a 
second order differential operator with constant coefficients and with the non-Hermitian 
Hamiltonian H first studied by Schwartz [12]. Then we give an example where 77 is a 
second order differential operator with variable coefficient. 



7.1. r] is a second order differential operator with constant coefficients 

Let us fix real constants d < and b (see (54)). In the simplest case potential V (x) 
may be a nonnegative constant Vo(x) = d 2 so that the solution of the spectral problem 
(31), (32) reads 

^\x) = )J^sm(kx), \{k) = k 2 + d 2 , k>0. 

Since H = —d 2 /dx 2 + d 2 , from equation (25) we find the absolute value of the 
transformation function 

p(x) = e dx . 

Next we solve equation (48) for u>(x) with 7 as given in (55) 

ui(x) = x = bx , b = . 

K ' 2d ' 2d 

(integration constant is without importance here) thus reconstructing the transforma- 
tion function (15) 

u(x) = e {d+ib)x . 

According to (20) function W(x) is constant W(x) = d. Therefore potential V(x) (23) 
coincides with V (a;) = d 2 . 

Transformation operator L p (27) is given by 

L p = —d/dx + d 

with the help of which we find eigenfunctions of H (18) 

* fc (x) = L p ^f(x) = ^-(k 2 + d 2 )- 1/2 [dsm(kx)-kcos(kx)}. 
Using unitary operator 

U = e~ lbx 
we obtain eigenfunctions 

ty k (x) =U^ k (x) = (k 2 + d 2 y 1/2 e- lbx [dsin(kx) - kcos(kx)} . 

of operator 

v = LL^=(~ + d-ib)(4- + d + ib)=~-2ib4- + ^ + dP. 
\ dx J \dx J dx A dx 
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Using equations (43) and (54) yields corresponding non- Hermit ian Hamiltonian 

H=-*-. 

dx 2 

Its Hermitian SUSY partner follows from (50), (49) and (55) 

h = -—. 

dx 2 

Note that although both H and h are defined by the same differential expression, 
corresponding spectral problems are different. For the Hamiltonian H we impose the 
boundary condition 

ip' k (0) + (d + ib)ip k (0) = (68) 

whereas for the Hamiltonian ho the boundary condition reads V'fc(O) = 0- Another 
remark we would like to make is that the spectral problem for H was first studied by 
Schwartz [12] as one of the simplest problems where the spectral singularity occurs in 
the continuous spectrum of H if d = 0. 

7.2. i] is a second order differential operator with variable coefficients 

According to Section 5, one can start with any scattering Hamiltonian ho with a real- 
valued potential Vq(x). The main element of the whole construction is the transformation 
operator L* which intertwines ho and a non-Hermitian Hamiltonian H (51). 
Let us choose 

2a 2 

v (x) = -o- -, a>0, c>0. 

cosh (ax + c) 

For c > the Hamiltonian h = —d 2 /dx 2 + vq(x) has no bound states in L 2 (0, oo). A 
solution of the Schrodinger equation with the Hamiltonian h having the asymptotic 
behavior as given in (54) reads 

e (d+ib)x 

u ( x ) = ; r [atanh(ax + c) — d — ib] . 

a — d — ib 

Using (13) and this transformation function, one finds function 

2a 2 

w{x) — d + ib-\ 27 r 

a sinh(2ax + 2c) — 2(d + ib) cosh (ax + c) 

which according to (10) and (53) defines both r\ operator and the non-Hermitian 
Hamiltonian H respectively. 

To illustrate selfadjointness of rj on Dh (4) we find operator H (18) unitary 
equivalent to rj. For that using (20), we first calculate W(x) 

. a 2 sech (ax + c)[ a tanh(aa; + c) — d] 

b 2 + d 2 + a tanh(ax + c) [ a tanh(ax + c) — 2d ] 
Then using (19), we find the potential 

A n 2( n 2_J2\ T2a 4 b 2 

V (x) = d 2 + d) - (69) 

W(x) W 2 (x) 
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where 

W(x) = b 2 + d 2 - a 2 + (a 2 + b 2 + d 2 ) cosh(2ax + 2c) - 2adsmh(2ax + 2c) . 

It is not difficult to see that V(x) — > d 2 as x — > oo by an exponential rule ~ exp(— 2ax) 
and, hence, as expected, operator H is scattering and positive definite. 

Another important point is that V(x) is continuous and bounded below. In 
this case, according to a known result [14], Operator H initially defined on a set of 
finite and twice continuously differentiable functions y(x) satisfying boundary condition 
y'(0) cos0 + |/(0) sin0 = has a closed selfadjoint extension H. Its domain where 
H = is described by the following properties [14]. If z(x) G -Djj then 

(a) z(x) G L 2 (0, oo), 

(b) z(x) is continuous and has an absolutely continuous derivative in any finite interval 
belonging (0, oo), 

(c) -z"{x) + V(x)z(x) G L 2 (0, oo), 

(d) z'(0) cos + 2(0) sin = 0, G M . 

Evidently this domain differs from that described in (22) only by the condition (b) which 
we everywhere skipped for simplicity. Going back from H to its unitary equivalent 
operator 77, we conclude that r\ = on D v = Dh (4). 

8. Conclusion 

In this paper basing on the factorization property of the metric i] operator [4], we 
proposed a special differential form for this operator. It happens that such a metric 
operator is unitary equivalent to the usual one dimensional Hamiltonian and therefore 
it has a natural supersymmetric structure. Once the metric operator is fixed, one 
can reconstruct the corresponding non- Hermit ian Hamiltonian. This opens a way for 
starting not from a non-Hermitian Hamiltonian and looking for the corresponding metric 
operator, but for going in the opposite direction, i.e., for starting with a given metric 
operator and presenting the corresponding non-Hermitian Hamiltonian. An advantage 
of this approach is the possibility to use the technique of SUSY QM for studying 
properties of the metric operator. In particular, as a possible line of future investigation, 
we are planning to use shape invariant metric operators and study how this property 
is reflected by corresponding non-Hermitian Hamiltonians. We pointed out that such 
metric operators are unbounded. Therefore some vectors from the initial Hilbert space 
are lost. They cannot be mapped into the new Hilbert space where the non-Hermitian 
Hamiltonian becomes Hermitian. 

We illustrated our approach by two examples. In the first example the metric 
operator is a second order differential operator with constant coefficients and w'(x) = 0. 
It happens that in this case the corresponding non-Hermitian Hamiltonian coincides 
with an operator proposed by Schwartz [12] (see also [13]) for illustrating the existence 
of the spectral singularity in the continuous part of the spectrum of this operator. In 
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the second example we considered a differential operator with variable coefficients as 
the metric operator. Using operator H (18) unitary equivalent to rj, we have shown that 
i] is selfadjoint on Dh (4). 
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